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OO ' We demonstrate that the classical dimer model defined on a toroidal hexagonal lattice 

. acquires holonomy phases in the thermodynamic limit. When all activities are equal the 

lattice sizes must be considered mod 6 in which case the finite size corrections to the bulk 
partition function correspond to a massless Dirac Fermion in the presence of a flat connection 
^ ' with nontrivial holonomy. For general bond activities we flnd that the phase transition in this 

model is a topological one, where the torus degenerates and its modular parameter becomes 
real at the critical temperature. We argue that these features are generic to bipartite dimer 
. models and we present a more general lattice whose continuum partition function is that of 

a massive Dirac Fermion. 



The classical statistical mechanics of lattice dimer models that interact only via hard-core 
exclusion are among the simplest exactly solvable two dimensional models. The subject has a long 
history and attracted the interest of both mathematicians [l| and physicists [2]. Recently, there 
, ■ have been several surprising advances and these models have reemerged in new and diverse areas 

ranging from quantum dimer models [1| to the melting corner /quiver gauge theory circle of ideas 
in string theory [4,]- Furthermore, the dimer model on the hexagonal lattice is closely related to 
^ ■ the physics of graphene [H] • 

. Dimer models divide into two classes: Those defined on bipartite and non-bipartite lattices — a 

lattice is bipartite if its sites can be coloured black and white so that adjacent sites are always 
of opposite colour. The most general model then assigns a positive number, called an activity, to 
each bond. 
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OO , In this note we focus on the simplest case of dimers on the hexagonal lattice, with three distinct 

activities a, b and c. We present two principal observations: 



• Surprisingly when all activities are set to one the model has a mod 6 structure and the limiting 
partition function is that of a free massless Dirac Fermion coupled to flat connections. 

• As the activities are varied there is a phase transition Q]- We show that at this transition 
the geometry collapses in a topological phase transition. 

The model is critical when none of the activities is larger than the sum of the other two. When 
considered on a torus, the finite size corrections are that of a free massless Dirac Fermion coupled 
to a flat connection. The angle between the periods of the torus goes smoothly to zero so that 
the torus collapses to a degenerate one dimensional structure at the transition. This geometrical 
collapse is surprisingly similar to that found recently in a very different context 0]. 

It has been found [H, 0] that the phase diagram of a bipartite dimer model is described by an 
amoeba obtained from the zero locus of a certain spectral polynomial, p{z,w), that arises in the 
solution of the model on a torus as described below. The amoeba constitutes the parameter domain 
where the model is critical. In general it has compact and bounding phase boundaries called ovals. 

Our results imply that the general situation for bipartite dimer models is therefore as follows: 
On the amoeba the partition function has finite size corrections corresponding to that of a single 
massless continuum Dirac Fermion on a torus, coupled to a flat connection. The dimer activities 
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and periodic repetitions of the fundamental cell determine the torus modular parameter and the 
holonomies of the flat connection. As the boundaries of the compact ovals are crossed there is a 
phase transition in which the Fermion acquires a Dirac mass [8]. 

The torus modular parameter, r, gradually changes as one traverses the amoeba until finally as 
the outer boundary is crossed there is another phase transition where the continuum geometry itself 
collapses. This transition is continuous — there is no latent heat — but asymmetric. As the transition 
is approached from the high temperature side the specific heat diverges with critical exponent 
a = 1/2, but there is no divergence as the transition is approached from the low temperature side. 

This transition is very similar to that discussed recently in a three matrix model 0]. There, as 
the critical temperature is approached from below, the spherical background geometry evaporates 
with a diverging specific heat and exponent a = 1/2; yet there is no divergence if the transition is 
approached from the high temperature non-geometrical side. 

We begin by constructing the Kasteleyn matrix for the N x M covering of a toroidal hexagonal 
lattice with fundamental tile given in Fig. [TJ We then review the properties of the bulk partition 
function and discuss the finite size corrections. The note ends with some speculations on the 
six-vertex model. 

If activities a = e"^"^" , b = e"^'^'' and c = e~^'''' are assigned to the bonds (c.f. Fig [T]) then the 
partition function of the model is given by 

Z{N, M, a, b, c) : 



E 



l^^b^^C 



(1) 



coverings 



where Ni is the number of active bonds of type i and Na + Nf^ + Nc = NM, since the lattice must 
be completely covered. When the activities are set to one, eqn. ([T]) counts the number of lozenge 
tilings of the dual triangular lattice. 





FIG. 1: The basic tile for the hexagonal lattice, showing the activities and our choice of Kasteleyn orientation and 
its amoeba. 



By assigning signs judiciously to the adjacency matrix, Kasteleyn j^, 10], Fisher and Temperley 
111 . Il2l ] observed that, one can convert ([T]) into a Pfaffian of, what is now called, a Kasteleyn 
matrix. 



The "clockwise odd rule" 13| is used: Arrows are placed on the links in such a way that 
by assigning +1 when following an arrow and —1 when opposing one, the product of the signs 
associated with any fundamental plaquette is —1 when the plaquette is circulated in a counter 
clockwise direction. 

The Kasteleyn matrix element K^b is then given by the activity assigned to the link from 
vertex A to vertex B times +1 if the arrow is from A to B but times —1 if it is from B to A. 

A Kasteleyn matrix, K, is therefore a signed weighted adjacency matrix for the lattice. On any 
simply connected planar domain, the modulus of the Pfaffian of K gives the partition function of 
the dimer model. 

On a toroidal lattice the partition function is given by a sum over the different discrete spin 
structures and therefore involves four terms. This Kasteleyn matrix can in fact be viewed as 
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a lattice Dirac operator [IJ], or more precisely a lattice version of CT> where C is the charge 
conjugation matrix and T> is the Dirac operator. In general a Kasteleyn matrix will describe 
many lattice Fermions with different masses and the theory of classical dimer models is essentially 
equivalent to that of two dimensional lattice Dirac operators. 

A consistent sign assignment for the hexagonal lattice is for all arrows to come out of the black 
vertices and into the white vertices. However, this assignment will not give the standard association 
of spin structures to the resulting Kasteleyn matrix. So we choose a slightly different convention 
with a and h bonds ingoing to the black vertex and the c bond outgoing. Then, for this hexagonal 
lattice, our choice of Kasteleyn matrix is 

K-i (2) 



where the matrix A becomes 

This matrix can be diagonalized using the eigenvectors z^w~^ and one finds 

AkMiM^'^'w-'^ = (c - a/z - bw)z^'w~'' (4) 
so the eigenvalues are given by the spectral polynomial 

p^'''^{z,w) = c- a/ z - hw . (5) 
Curiously, this same polynomial also determines the phase diagram in Kitaev's quantum honeycomb 



model [15 1 



For a more complicated tiling (c.f. Fig. [2|) the analogue of ^ will result in a Fourier transformed 
Kasteleyn matrix whose off diagonal block is of size d where d is the number of vertices of the same 
colour in the fundamental domain. The determinant of this matrix will again be a polynomial in 
z and w which will determine the phase diagram for the corresponding bipartite dimer model. On 
the amoeba of this polynomial K describes one massless Dirac Fermion and d — 1 massive ones, 
with the massive modes decoupling from the problem in the thermodynamic limit. 

In the thermodynamic limit the logarithm of the bulk partition function per dimer, W = 
is found to be 

^^-rSrS'"(.-»e-»-.e-*). (6) 

If one of the weights is larger than the sum of the other two, then W is the logarithm of that 
weight. There are three such regions, called frozen regions, where one has Vl^ = Ina, W = lnb and 
W = Inc |lo| . The region where none of the weights is larger than the sum of the other two is 
referred to as the amoeba of the spectral curve p^'^^{z,w) and is shown in Fig. [TJ On this curve 
there are exactly two pairs of angles (O, <I>) and (—0, at which p^'^^ = where, 

sin(e) = ^ , sin($) = ^ and (7) 
2r 2r 

abc 

r = ^ (8) 

Y^(a + b + c){-c + a + b){c - a + b){c + a - b) . 

Here, r is the radius of the circumcircle of the triangle with sides a, b and c and opposite angles 
and TT — — respectively. The radius r diverges, while and $ become zero or vr, on the 
boundaries of the amoeba. 
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A little effort sliows tliat W{a, b, c) can be expressed in the form 

Wia,b,c) = lnc + -ln(6/c) + ^{U2{-e'®) - /^aC-e"'®)) (9) 
vr Im c c 

wliere it was assumed tliat c>h> a. Using the Lobachevsky function ^ 

L{z) = - [ ln|2sin(t)|(it (10) 
Jo 

W can be written in the more symmetrical form 

W{a, b, c) = Palna + phlnb + Peine + ent{a, b, c) (11) 
where ent(a, b, c) is the entropy per dimer and is given by 

ent{a, b,c) = - {L{TTpa) + L{Trpb) + L{ttPc)) (12) 
vr 

with Pa + Pb + Pc = ^ and pa = ^/n and pb = Q/tt. 

The quantities pa, Pb and pc are the probabilities that the bonds with activity a, b and c 
respectively, belong to a randomly chosen dimer configuration 17l |. 

As the edge of the amoeba associated with say, c, is approached pc — > 1- Curiously ent{a, b, c) 
is the Bloch-Wigner dilogarithm D{z) where z can be chosen to be the ratio of any two activities 
(such that \z\ < 1) times the exponential of i times the angle between them. This latter quantity 
has iC-theoretic significance 18l |. 

The internal energy U = — is continuous throughout the phase diagram (see [i3) while the 

specific heat C = P"^^^ diverges at the transition Q]. 

Explicitly, for a = b = and c = 1, for /3 > ln2 we have W = U = C = 0, and for < ln2 
we have 

2 /-In 2 y 

W=- cos-i(-)dy, (13) 

J/3 2 

2 . e>^ B"^ e^ 
U = -cos-\—) and C = - , (14) 

Expanding around /3c = In 2, with In 2 — /? > we have 

4\/2 3 

W{(3)^^{\n2-0rK (15) 

In sum: The transition is continuous with no latent heat. The specific heat is zero in the low 
temperature frozen phase; there is a phase transition at /? = In 2, and the specific heat diverges 
with critical exponent a = ^ as the transition is approached from the high temperature side. 

On the torus corresponding to the N xM covering by the fundamental tile the partition function 
is given by 

(_^7VM(JVM-l)/2i ^2mil+u+v+2uv)p^^Qj^^^^ 
u,v=0 

1 

= ^Y1 e2"^(5+«+«+2™)DetA,« (16) 

u,v=0 



^ L{z) — I (7/2(22:), where Chix) is Claussen's function. 
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with the four terms corresponding to the four discrete spin structures on the torus. For the tihng 
in Fig. [T] 



N-l M-1 

27ri{n-\-u) 27vi{m-'rv) 



n-m r rr ZTTll n-HTJ, I Z7VZ[ m-f-V I 

Hp'^'^ie W,e^M^) 

m=0 

((c-6e^^%^)^-a^e2™) . (17) 



71=0 m=0 
AI-1 



n 



m=0 

The polynomial p^'^^ has two zeros: One at z = e*® , = e** the other at z = e~*® , = e~**. 

Consider, for simplicity, the case when the weights a, b, and c are such that the angles Q and $ 
are rational multiples of 27r, and in particular such that = — and $ = , for integral ng 
and rriQ. Then shifting n and m appropriately and expanding around the zero at (0, <1>) we have 

p"^"" -^^e-'^{n + u + T{m + v)) (18) 

where 

r = i^e^(e+*) . (19) 

Ma ^ ' 

In the scaling limit, A^, M — > oo, with the activities and ^ = ^ held fixed, r remains constant 
and becomes the modular parameter of the continuum limit torus. The second zero at (—0, — $) 
simply changes r to r. 

Then for large M and N we have 

r Z{N,M) 



1/2 

2 ^ 

^^,^!=0 



"J(0|r) 



?7(r) 



(20) 



Where Z£)iraciT,0,(l)) is the partition function for a Dirac Fermion propagating on the continuum 
torus with modular parameter r in the presence of a gauge potential with zero field strength, i.e. 
a flat connection, but with holonomies e^'^*^ and e'^'^^'f' round the cycles of the torus. 

The two zeros of p^'^^{z,w) guarantee that we are dealing with a continuum limit for a Dirac 
Fermion, and in the current tiling the Fermion doubling phenomenon is minimized since there are 
only two vertices in the fundamental tile. A rectangular toroidal domain arises when = + 6^. 
which corresponds to tq = and ti = Cf- 

For infinite temperature, where a = b = c = 1 one can see that, if N and M are positive integral 
multiples of 6 then Det^o,o = in (I17p but nonzero otherwise. A similar structure can be seen for 
other values of u and v. A careful study of the limiting form of (|17p therefore requires that the 
limit be taken with the mod 6 periodicity in N and M made explicit. Setting N = p mod 6 and 
M = q mod 6 we obtain ([20]) with the modular parameter r = , ^ = ^ — | and </> = ^ + f- 

As the temperature is decreased the angle between the periods of the torus, given by the phase 
of T, reduces continuously reaching zero when the boundary of the amoeba is reached. For the 
special case a/c = b/c = the modular parameter takes the form 

r~e(l-4(ln2-/3)+i2\/2Vln2-/?) +... . (21) 



We see that, at the critical temperature Pc = In 2, where ri = while tq = ^, the torus degenerates. 
It effectively reduces to a line. This is a singular limit of the geometry and is quite distinct from 
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the cyUnder hmit where for a torus with fixed angle both tq and ri go to zero. This phenomenon, 
of a singular limit as the boundary of the amoeba is approached, is generic. 

Note: For /3 = we have t{^^) = ^ + i^ = ~^fjWj- With zero holonomies, the finite size 
corrections are modular invariant, and so unchanged under the replacement r — > — ^, however, the 
bulk term is different at these two temperatures. 

Note also, that special care is again needed for and $ rational multiples of 2tt. E.g. when 
a = b and = with p and q relatively prime, so that /? = ln(2cos(0)), then the lattice sizes 
M and N must be considered modulo q and one encounters a rich structure of flat connections 
contributing to the finite size corrections. 

The dimer model is also closely related to the six vertex model and in fact on the free-Fermion 
line the six vertex model reduces to a dimer model. In the disordered phase the zero field six- 
vertex model is critical [l^ . As the parameters of the six- vertex model change there is a phase 
transition, which generically has a latent heat, to a ferroelectric frozen phase [20]. But when the 
phase boundary is crossed from the disordered phase the specific heat diverges, again with exponent 
a = 1/2. 

We conjecture that, for the six-vertex model on a torus, as the ferroelectric phase is approached 
from the disordered phase, the modular parameter, r, decreases to a critical but non-zero value, 
Tc, and a topological transition occurs. The special case of the dimer model is consistent with this 
and has Tc = consistent with the free-Fermion line. 

It was found in ^2lj] (see also [23]) that the six- vertex model with domain wall boundary con- 
ditions could be rewritten as a matrix model. The spectrum of eigenvalues of this matrix model 
changes from continuous to discrete as one crosses the critical line between the disordered and 
ferroelectric phases. Though this matrix realization of the six vertex model has not been extended 
to toroidal boundary conditions it gives some idea of what is happening at a microscopic level as 
the transition is crossed. 

The situation is again reminiscent of what happens in the fuzzy sphere model 0]. However, 
the high temperature phase is geometrical with a continuous spectrum in the six-vertex and dimer 
models, while the geometrical phase in [6,] has discrete spectrum and occurs at low temperatures. 
The common feature is the asymmetruc structure and characteristic divergence of the specific heat 
in all of these models. We expect that this is generic to such topological transitions. 

Finally, we observe that the full generality of features of the bipartite class of dimer models 
can be accessed in the model whose fundamental tile amd amoeba are shown in Fig[2j This figure 
shows a tiling of the torus whose phase diagram or amoeba contains a single compact oval. 





FIG. 2: Four copies of a Kasteleyn oriented basic tile for the tunable single oval model and its amoeba with 
A± ^ B± = C± = 1 and D = 10 in Eq. 



The corresponding spectral polynomial is in general 



p 



Oval 



[z, w) = D A^z 

z w 



w 



z 



(22) 



where the seven parameters, A±, B± etc., are all positive and determined by the twelve link 
activities and can be parametrized as D = 2Acos\i{tx) + 2B cos\i{ty) + 2C cosh.{txy), A± = Ae^*^, 
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B± = Se^*™, and C± = Ce^^". As the oval is crossed the continuum partition function describes a 
massive Dirac Fermion and both the Dirac mass and modular parameter are tunable by adjusting 
the activities 
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